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ABSTRACT 



An investigation was made into the types of fixing error that 
would result when two passive directional sensors in one velocity 
medium was obtaining its fix on a target in another velocity medium. 
For typical velocity differences that could exist in the oceans off 
the U. S., this error was investigated for various values of sensor 
spacing and target/sensor relationships to a velocity discontinuity. 
Tabular and graphical results were obtained which revealed that 
errors could result under certain conditions and these results were 
predictable and could operationally be avoided. 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5101 



TABLE OF CONTENTS 



Section Page 



1 . 


Introduction 


9 


2. 


Problem Modeling 


13 


3. 


Mathematical Solution 


18 


4. 


Example 


22 


5. 


Results 


23 


6. 


Conclusions 


31 


7. 


Bibliography 


50 


8. 


Appendix I 


51 



3 



v ,■ i ij xo i > y.' f 

' ' ».nnv q • 

-t . 'J,r -• i 



List of Tables 



Table 



1. 


Relative 


Error 


< v i 


= 


1500 


2. 


Relative 


Error 


(Vi 


= 


1525 


3. 


Relative 


Error 


(Vi 


- 


1525 


4. 


Relative 


Error 


(Vi 


s 


1500 


5. 


Relative 


Error 


(V X 


= 


1550 


6. 


Absolute 


Error 


(Vi 




1500 








(Vi 


= 


1525 





Page 


V 2 = 1525) 


39 


V 2 = 1500) 


41 


V 2 = 1475) 


43 


V 2 = 1550) 


45 


V 2 = 1500) 


47 


V 2 = 1525) 


49 


V 2 = 1500) 





5 



Figure 

1 . 



LIST OF ILLUSTRATIONS 



Page 

11 



Formation and Course of the Gulf Stream 



2a. 


SST 8-12 January 1964 




12 


2b. 


Approximations of Sea Surface Temperatures 


13 


3. 


Horizontal Ray Bending Due to Different 
Velocity Media 


15 


4. 


Snell's Law (V L > V 2 ) 




16 


5. 


Illustration of Critical Angle 




17 


6. 


Changes to Problem if Sound Source is 
Quadrant Four 


in 


21 


7. 


Relative Error in Miles for A/B Ratio 
1, V L = 1500m/s, V 2 * 1525m/ s 


of 


25 


8. 


Relative Error in Miles for A/B Ratio 
1, V L = 15 25m/ s, V 2 = 1500m/s 


of 


26 


9. 


Relative Error in Miles for A/B Ratio 
5/1, V L = 1500m/ s, V 2 = 1525m/s 


of 


27 


10. 


Relative Error in Miles for A/B Ratio 
5/1, V 1 = 15 25m/ s, V 2 = 1500m/s 


of 


28 


11. 


Relative Error in Miles for A/B Ratio 
1/5, V L = 1500m/s, V 2 * 1525m/s 


of 


29 


12. 


Relative Error in Miles for A/B Ratio 
1/5, V L = 1525m/ s, V 2 = 1500m/s 


of 


30 


13. 


Absolute Error in Miles for A/B Ratio 
1, V L = 1500m/ s , V 2 = 15 25m/ s 


of 


32 


14. 


Absolute Error in Miles for A/B Ratio 
1, V L = 1525m/s, V 2 = 1500m/ s 


of 


33 


15. 


Absolute Error in Miles for A/B Ratio 
5/1, V 1 =■ 1500m/ s , V 2 = 1525m/s 


of 


34 



7 



\ 



F igure 
16. 


Absolute 


Error in Miles 




5/1, 


Vl 


= 1525m/ s, V 2 


17. 


Absolute 


Error in Miles 




1/5, 


V 1 


= 1500m/ s, V 2 


18. 


Absolute 


Error in Miles 




1/5, 


V 1 


= 15 25m/ s, V 2 


19. 


Relative 


Error in Miles 




i. 


V 1 - 


1500m/s, V 2 = 







Page 


for A/B Ratio 
= 1500m/ s 


of 


35 


for A/B Ratio 
= 1525m/s 


of 


36 


for A/B Ratio 
= 1500m/ s 


of 


37 


for A/B Ratio 
1550m/ s 


of 


38 



8 



1 . 



Int roduction. 



Though fixing by means of cross bearings from two or more 
directional acoustical sensors is an established concept in Naval 
Warfare, there has been little investigation done as to the effects 
of fixing error due to horizontal water gradients. With the advent 
of the SSQ-36 sonobuoy a tool is now available to the ASW type air* 
craft that allows measurement of the temperature of the water, 

Sound waves travel in straight lines only in a medium in which 
the speed is everywhere constant. The speed of sound in sea water 
depends on salinity, pressure, and temperature. The problem of 
salinity and pressure can be neglected when discussing horizontal 
gradients in open ocean environments since both will be almost 
constant throughout the media. With decreasing temperature, the 
water becomes more dense, and the velocity of the sea water de- 
creases. A useful equation for relating these variables is 

V ft /sec = 4502 + 9.184T - 0.0289T 2 + 0.0182Y (1-1) 

where Y * depth in positive feet 
and T * temperature. 

Using (1-1) gives a velocity of about 4900 ft/sec or 1500 m/sec 
for temperatures of 50°F. Similarity, a temperature of 70°P and 90°F 
gives about 1525 m/sec (5000 ft/sec) and 1550 m/sec (5100 ft/sec) 
respectively. 

An operational problem was envisioned in which a submarine was 
trying to penetrate a barrier that had been established either in 
or outside of the Gulf Stream. The submarine was then outside of, 
or in the Gulf Stream. This would imply a very distinguishable 
temperature gradient existed between the buoys and the submarine 
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which would imply that a distinct horizontal velocity gradient existed. 

In the ocean are various currents. Many of these currents are of 
a different temperature than the water they are flowing through. At 

this point a digression will be made to give a short explanation of 

c v. v. > 

currents to the reader. Although it is not absolutely necessary for 

j T * . ^ l i i .. 

understanding the problem of horizontal gradients and its solution, 

t- ' i l • * s ,/. 

it will allow a better understanding of the mechanisms involved. 

The two primary causes of currents are the exchange of heat 
with ’the atmosphere and energy supplied by winds. Even in 
the absence of winds, a thermal circulation would be set up 
because heating in lower latitudes would warm the surface 
water, causing them to expand slightly, whereas cooling in 
high latitudes would cause them to shrink. The waters in 
low latitudes would therefore stand a little higher and there 
would be a flow away from the equator towards higher lati- 
tudes. It is believed that the winds have a much more im- 
portant effectl Not only do they exert a drag on the sur- 
face of the sea, setting the surface layers in motion but 
they also will bring about a redistribution of mass 
establishing slopes. /J 2 J 

However, the intense currents such as the Gulf Stream have been 

> 1 J *• ■ : :' J i i' l « 

shown not to be related to local winds. It results from the varia- 

* * ( * . “J * i . 

tion of Coriolis Force with latitude. The winds only add energy. 

'-'.vs vr. 

As an example of the size of a major current, the Gulf Stream starts 
out in Florida about 95 nautical miles wide and nearly 2 miles deep 

. t ' . * 

with a speed of 3.5 knots. Off Cape Hatteras it has been reduced 

to 1 knot, and becomes wider, more shallow and more interrupted as 

• 0 1 l 

it moves out into the Atlantic. The average temperature of the 
surface water is 80°F as far north as Cape Hatteras and the eastern 
side is warmer than the western. Figure 1 shows the Gulf Stream 
and its different component currents. Figure 2a shows the intricate 
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isotemps of sea surface temperatures for a four day period. To 
simplify the problem a first order approximation is made as 
illustrated in Figure 2b. The temperatures of the Gulf Stream 
are broken down into blocks or strips. 



\ 




Formation and Course of the Gulf Stream 
Figure 1 
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Boston 




SST 8-12 January 1964 
Figure 2a 
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Approximations of Sea Surface Temperatures 
Figure 2b 



2. Problem Modeling, 

When structuring the problem the following assumptions were made: 

(1) An initial detection was already given; thus a detection 
range limitation was not considered. It made no difference 
if detection was the result of bottom bounce, channeling, 
direct beam, or surface reflection. This allows vertical 
gradients to be neglected since they have almost no effect 
on the problem when a detection is given. The difference 
between horizontal gradients at various depths does not 
change appreciably. 

(2) A directional acoustic sensor was considered to have no 
bearing error, i,e., finite beamwidth. Adding this problem 
would be the next order of sophistication. 
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(3) Although both horizontal continuous and discontinuous grad- 
ients occur in the ocean, the temperature gradients were 
considered sharp and distinct, i.e., discontinuous. 

Figure 3 exhibits a pictorial display of the problem. Two 
sonobuoys are placed 2 a apart from each other and perpendicular to 
a distinct temperature gradient which is a distance b away. The 
origin ( 0 , 0 ) of the problem is the midpoint between the two sono- 
buoys. Consider the sound source actually at point (X 3 ,Y 3 ) . For 
buoy 1 the sound ray path starts with an angle 02 to the normal 
and when it reaches the gradient, the ray bends toward the normal 
and becomes 0j_. Similarily for buoy 2 the ray path starts with an 
angle 04 to the normal and then is bent to become 03 . 




Figure 3 
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When a sound source must pass through a velocity differential 
the acoustical Snell's Law can be applied to understand what happens. 
Snell's Law says 

Sin 6j = Vi ( 2 - 

Sin 0 t V 2 

where 0^ = angle of incident ray 

0 t = angle of transmitted ray 
V = velocity of medium. 

If is greater than V 2 Figure 4 is applicable. This implies that 
the angle is bent toward the normal or < 0^. 




Snell's Law (V^> V 2 ) 

Figure 4 

If V 2 is greater than V^, (i.e., < V 2 ) then Figure 5 is applicable 

and we have the unusual phenomenon which is described as a critical 



angle 0 C occurring when 
Sin 0 C = 

V2 
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(2-2) 




Illustration of Critical Angles 
Figure 5 

Here the angle is bent away from the normal or 0 t ^ 0 j. until 0 C 
is reached at that point and thereafter the incident ray is bent 90 ° 
(critical angle) and no sound (theoretically) is transmitted intQ the 
second velocity medium. In reality some sound is transmitted. How- 
ever, it is so small that it will be neglected. 

If it is assumed that the reader is at the sensor, the rays 
enter with angles «i and 83 (Figure 3 ) . Projecting their path back- 
ward the sound source appears to be at the point (X4, Y4) . Here the 
error involved is arbitrarily selected and is the vector from (X^.Y^) 
to (X3,Y3). The error angle 85 is taken from the apparent extention 
of the ray with angle 83 and the error vector. Finally if the case 
arises where X3 ^ b, then (X3,Y3) « (X^,Y^) since ■ V2 and no 



error exists. 
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3. Mathematical Solution.' 



Now that the actual and apparent sound paths have been shown, it 
is necessary to derive the accompanying equations in order to solve 
the problem. For reference see Figure 3. 



Let (X 2 ,Y 2 ) and (X^,Y;l) be the points of intersection at b of 
ray path 1 and ray path 2 respectively. Equation of line from 
(0,-a) to (X 2 ,Y 2 ) is , - 





Y = X Tan 9i -a 


(3-1) 


or since X 2 = b, 


Y 2 = b Tan 9^ -a 


(3-la) 


Equation of line 


from (X 2 ,Y 2 ) to (X 3 ,Y 3 ) is of the form 






Y-Y 2 = (X-X 2 ) Tan 9 2 


(3-2) 


Combininb (3 -la) 


and (3-2) 






Y-b Tan 9i + a = (X-b) Tan 92 


(3-3) 


Previously it was 


shown that Snell’s Law says,,.. 






0 2 = arcs in (^2 s in 9 i) 

Vl 


(3-4) 


Combining (3-3) and (3-4) gives the equation of line A 

4 , L’> ( . 




Y-b Tan 9 L ' 4 


V 2 

■ a = (X-b) Tan (ar'csin ( — sin 9^)) 

V 1 


(3-5) 


✓ « 

Similarily the equation of line from (0,+a) to (X 25 Y 2 ) is of 




the form 


Y = X Tan 93 4- a 


(3-6) 


or since X 2 = b, 


Y 2 = b Tan 93 4- a 


(3 -6 a) 


and the equation 


from (X]^) to (X 3 ,Y 3 ) is 






Y-b Tan 9 3 -a = (X-b) Tan 9 4 


(3-7) 
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(3-8) 



Combining (3-6a) and (3-7) gives the equation of line A' 

Y-b Tan 83 -a * (X-b) Tan (arcsin ( ^2 fl £ n 

The apparent intersection (X^jY^) for 83 > 0 is the inter- 
section of (3-1) and (3-6) or 

X 4 Tan 8 ^ -a * X 4 Tan 83 + a (3-9) 

Solving for X 4 , 



XA = 2 a 

Tan 0^ - Tan 83 


83 > 6 


ombining (3-1) and (3-6) 




Y 4 Tan 83 + a Tan 83 * Y 4 Tan 8 j - 


a Tan 9^ 


y _ a (Tan 83 + Tan 0i) » 


e 3 ^ 0 


Tan 0^ - Tan 83 



(3-10) 

(3-11) 

(3-12) 



The real position (X 3 ,Y 3 ) for 83 > 0 is given by the inter- 
section of A and A'. Rewriting equation (3r5) and (3-8), 



Y * b Tan 8^ -a + X Tan (arcsin(X2. sin 8^)) -b Tan (arcsin(X2 sin 0j)) 

V 1 (3-13) 



Y = b Tan 83 + a + X 



Tan (arcsin (^2. sin 83 )) -b Tani (arcsin (^2 sin O 3 )) 

1 1 (3-14) 



To facilitate writing equations, the following definitions are 
added , 



Z^ ■ arcsin (X 2 . sin ( 0 j)) 



Z 2 * arcsin (^2 sin ( 83 )) 
V 1 

T l * tan(Zj) 

T 2 * tan(Z2) 
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Combining (3-13) and (3-14) and solving 
x 3 (T^-T 2 ) = b (Tan 9 3 - Tan 9 L ) + 2a + b (T L -T 2 ) 

_ b(Tan 9o - Tan 9i + Ti-To) + 2a. 

*3 “ > 

Ti-T 2 



Again combining (3-13) and (3-14) and solving 

X _ y 3 + a " b Tan ®1 + b T^ _ Y 3 - a - b Tan 93 + b T 2 
Ti T 2 



y 3 = b(T 2 Tan - T L Tan 9 3 ) - a (T t + T 2 ) 0 Q 

T 2-Ti ’ 3 



For the case 9 3 0, A is the same, i.e,, 

»_ 1 * • * 

Y 2 = b Tan 9-^ - a 

Y - b Tan 9 L + a = (X-b) T x 
However A 1 becomes 

Y 0 = - b Tan 9q + a 

i , ,u l $ 

l • • M 

Combining (3-19) and (3-7) gives the equation of A 1 1 

{• i • * .» 

Y + b Tan 9 3 - a = (b-X) T 2 

Using (3-1) and (3-19) gives the apparent intersection (X^,Y^) 
for 9 3 0 . 

X^ Tan 9]_ - a =-X Tan 9 3 + a 

X A = — , 9o< 0 

4 Tan 9 l + Tan 9 3 J 

Y 4 + a = a - Y 4 
Tan 9-^ Tan 9g 



(3-15) 

(3-16) 

(3-17) 

(3-la) 

(3-18) 

(3-19) 

(3-20) 

(3-21) 

(3-22) 

(3-23) 
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(3-24) 



Y = a (Tan 0 L - Tan 03 ) 


0 3 < 0 % 


Tan 0^ + Tan 0 3 





Solving now for the real intersection (X^.Y*) using (3-18) and (3-20), 



= b Tan 0^ - a + X^T^ - bT^ = -b Tan 63 + a + bTg - X 3 T 2 



X, = b ( T 2 " T l " Tan ®3 ' Tan ®l) + 2a . 0„ < 0 




Ti+T 2 3 


= a + b (T^ - Tan 9^) _ a-fb(T 2 * Tan 0^) • Y r 


) 


T 1 t 2 


Y, = a(T l - T 2 ) + b ( t 2 Tan ®1 ' T 1 Tan 0 3> , 0,< 0. 


> 


T 1 + T 2 





(3-25) 

(3-26) 



(3-27) 

(3-28) 



These equations are for the first quadrant. If the sound source 
is in the fourth quadrant make the following changes to Figure 3 as 
shown below in Figure 6 . Also change and 0^ to correspond. 




After change 




Changes to Problem if Sound Source is in Quadrant Four 

Figure 6 
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Now that all the necessary equations have been developed, it is 



now possible to find the relative percent error as a function of 
a, b, 9-p and Qg. This is done by taking the ratio of radial error 
or error vector to the actual distance of the sound source from the 
origin (0,0)- Relative error (RE) is 



Appendix I contains the computer program from which solutions 
for relative error were obtained. 

4. Example. 

To demonstrate the fo rementioned, the following example is worked. 
If two sonobuoys are placed on the inner edge of the Gulf Stream and 
the submarine is outside the Stream so that = 1525 m/sec and V 2 = 
1500 m/sec, a = 1 NM = b and 9-^ = 30°, 93 = 15°. Using equations 
(3-10) and (3-12) it is possible to solve for and respectively, 
where 



RE 




(3-29) 




X 



4 ~ -57735 - .26795 







Y = Cl) (.57735 + .26795) 
4 .57735 - .26795 



•8453 = 2.73 

.3094 



and the apparent position is (6.46,2.73). 
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Using equation (3-15) and (3-16) it is possible to solve for 
X 3 and Y 3 respectively, where 
Z x = .51416 
H 2 = .25741 
T 1 = .56483 
T 2 = .26325 

( 1 ) (.26795 - .57735 + .56483 - .26325) + (2) (1) „ 2 - .00782 

X 3 = .56483 - .26325 .30158 

= 6.61 

(1) (.26325) (.57735) - (.56483) (.26795) - (1) (.56483 + .26325) 

Y 3 = - .30158 

( .15199 - .15135) - .82808 _ 2 . 

- .30158 

and the actual position is (6.61, 2.74). 

To figure relative error, the relationship 

RE = 1 ( 6.61 - 6 . 46 ) 2 + (2.74 - 2.73 ) 2 = 

](6.61) 2 + ( 2 . 74) 2 

is used, 

5. Results, 

Basic results are graphically illustrated in Figures 7 through 
12 and tabulated in tables 1 through 5 for varying 0 ^ and 83 , 
different and V 2 and varying A/B ratios which show relative error. 
Percent relative error is plotted versus 83 , for given 9^, velocities 
and a/b ratios. In general it can be observed that: 



23 



(a) a/b is invariant; that is, only the ratio of a/b effects 
the results, not the value of a or b alone. For example, if a is 

1 mile and b is 2 miles or if a is 50,000 yards and b is 100,000 
yards the relative error is the same. 

(b) it was noted that there was not a large difference in 
error whether V-^ = 1500 and V 2 = 1525 or if V-^ = 1525 and V 2 = 1500. 
However, there was a consistently smaller relative error for V-^ > V 2 
compated to ^2^ ^ 1 * 

(c) for a ratio of a/b = 1 , it was noticed that the error 
was independent of 0^. This is not true for a/b either less than or 
greater than 1. For 63 — 0° the relative error was about 1.5% which 
gradually increased to 3% for 83 = 30°. Over 30° error increased 
fairly rapidly up to about 10% for 83 = 55°. At around 70° the angle 
became critical for V-^ = 1525 and V 2 = 1500. 

(d) for a/b 1 , for example a/b = 2 , sonobuoys spaced 20 
miles apart and the buoys placed 5 miles from the gradient, the error 

increased as 0-^ increased holding 83 constant. The error also in- 

» ^ 

creased as 83 increased. As the ratio of a/b increased the error 
increased. It should be noted that around 83 = 25° the error started 
increasing fairly rapidly. The error at its extreme reached over 12%. 
Average error appeared to be around 3%. 

(e) for a/b < 1, the reverse of a/b 1 was observed. As 
..." ! . t ^ 

0 -^ increased, the error became smaller except at the maximum extremes 
* * >• . . , o > 

for 0j and 83 . That is, unless 0^ - 83 < 10° as 0^ increases, error 

is smaller. A major reason for this is that as a/b-^0, most of the 

apparent positions fall in V^, so that error is 0 . 
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Relative Error in Miles 



Relative Error in Miles for A/B Ratio of 
1, V L =* 15 00m/ s, V 2 - 1525m/s. 




©2 in Degrees 
Figure 7 
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Relative Error in Miles 



Relative Error in Miles fox A/B Ratio of 
1, V 1 = 1525m/ s , V 2 = 1500m/s. 



5 



4 



3 



2 



1 



0 




-40 -20 0 20 40 60 



©2 in Degrees 
Figure 8 
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Relative Error in Miles 



Relative Error in Miles for A/B Ratio of 
5/1, Vj = 1500m/ s, V 2 - 1525m/s. 




«oor S 3QFlgure 9 

Oi siL'gi'i 



r<£ 
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Relative Error in Miles 



Relative Error in Miles for A/B Ratio of 
5/1, V L = 1525m/s, V 2 = 1500m/s. 




83 in Degrees 
Figure 10 



28 



Relative Error in Miles 



Relative Error in Miles for A/B Ratio of 
1/5, V L » 1500m/s , V 2 - 1525m/s. 




9^ in Degrees 
Figure 11 
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Relative Error in Miles 



Relative Error in Miles for A/B Ratio of 
1/5, V = 1525m/ s , V 2 = 1550m/s. 




9^ in Degrees 
Figure 12 
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Although relative error allowed an effective way to represent 
graphically the fundamental relationship, it is also important to 
look at the value of absolute error. This error is tabulated in 
Table 6 and in Figures 13 through 18. In general it can be observed 
that: 

(a) for a/b £ 1 there are no errors of more than 1 NM, 

(b) for a/b > 1 the absolute error is around .75 NM or 

less for 93 / 0 ^ 5 1/2. If 83/01 > 1/2 the error goes up to 5 NM. 

Figure 19 is inserted to illustrate what happens to relative error 
if the velocity is doubled. 

6. Conclusions. 

Based on theoretical results it can be concluded that: 

(a) Errors of 5 % or more never take place where angles less 

than 40 degrees normal to the temperature discontinuity (13° difference) 
occur. 

(b) Sonobuoy spacing can affect the realm of errors possible. 

In general, a small sonobuoy spacing will keep error to a minimum. 

(c) The farther sonobuoys are kept from a temperature discontinuity, 
the smaller the relative fixing errors will be, though the absolute 
error will not necessarily decrease. 

(d) If a reasonable model of the water's velocity structure is 
known, the true position of the target could be obtained. 
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Absolute Error in Miles 



Absolute Error in Miles for A/B Ratio of 
1, = 1500m/s, V 2 = 1525m/s. 




83 in Degrees 
Figure 13 
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Absolute Error in Miles 



Absolute Error in Miles for A/B Ratio of 
1, = 1525m/ s, V 2 = 1500m/s. 
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Absolute Error in Miles 






Absolute Error in Miles for A/B Ratio of 
5/1, V x = 1500m/ s , V 2 = 1525m/s. 




9^ in Degrees 
Figure 15 
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Absolute Error in Miles 



Absolute Error in Miles for A/B Ratio of 
5/1, » 1525m/s, V 2 - 1500m/s. 




-40 -20 0 20 40 60 

0^ in Degrees 



Figure 16 



Absolute Error in Miles 



Absolute Error in Miles for A/B Ratio of 
1/5, V l = 1500m/s, V 2 = 1525m/s. 




0^ in Degrees 
Figure 17 
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Absolute Error in Miles for A/B Ratio of 

1/5, Vi * 1525m/s, V 2 * 1500m/s. 




0^ in Degrees 
Figure 18 
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Relative Error in Miles for A/B of 
1, Vj_ “ 1500m/s, V 2 = 1550m/s. 




9^ in Degrees 
Figure 19 



38 



TABLE 1 

RELATIVE ERROR 
V t = 1500 V 2 - 1525 



9, =10° 0 3 




A/B=l 


A/B-l/2 


A/B=l/5 


A/B-l/10 


A/B-2/1 


A/B-5/1 


1 l.~ 


-40° 


1.24% 


0% 


0% 


0% 


1.89% 


2.28% 


2. 


-30 


1.28 


.50 


0 


0 


1.68 


1.92 


3. 


-20 


1.33 


.83 


0 


0 


1.58 


1.73 


4. 


-10 


1.41 


l.ll 


.20 


0 


1.57 


1.66 


5. 


- 5 


1.48 


1.25 


.57 


0 


1.59 


1.66 


6. 


0 


1.56 


1.41 


.95 


.20 


1.64 


1.68 


7. 


5 


1.66 


1.59 


1.35 


.96 


1.70 


1.73 



e,= 2 o° 8. 


-10 


1.33 


.83 


0 


0 


1.58 


1.73 


9. 


- 5 


1.41 


1.00 


0 


0 


1.62 


1.75 


10. 


0 


1.52 


1.18 


.16 


0 


1.69 


1.79 


11. 


5 


1.64 


1.38 


.58 


0 


1.78 


1.86 


12. 


10 


1.79 


1.61 


1.04 


.12 


1.89 


1.95 


13. 


15 


1.98 


1.88 


1.57 


1.07 


2.03 


2.06 



ei=30°14. 


- 5 


1.39 


.70 


0 


0 


1.75 


1.96 


15. 


0 


1.52 


.90 


0 


0 


1.83 


2.02 


16. 


5 


1.67 


1.12 


0 


0 


1.94 


2.11 


17. 


10 


1.84 


1.37 


0 


0 


2.07 


2.21 


18. 


15 


2.04 


1.66 


.54 


0 


2.23 


2,35 


19. 


20 


2.29 


2.01 


1.18 


0 


2.43 


2.52 


20. 


25 


2.60 


2.44 


1.98 


1.21 


2.68 


2.72 



0i =40°2l . 


0 


1.53 


..42 


0 


0 


2.09 


2.43 


22. 


10 


1.89 


.94 


0 


0 


2.37 


2.66 


23. 


20 


2.39 


1.64 


0 


0 


2.77 


3.00 


24. 


30 


3.14 


2.66 


1.23 


0 


3.39 


3.53 


25. 


35 


3.68 


3.39 


2.55 


1.18 


3.82 


3.90 



0i=5O°26. 


0 


1.41 


0 


0 


0 


2.47 


3.12 


27. 


10 


1.83 


0 


0 


0 


2.80 


3.39 


28. 


20 


2.39 


.69 


0 


0 


3.26 


3.79 


29. 


30 


3.22 


1.79 


0 


0 


3.96 


4.41 


30. 


40 


4.57 


3.56 


.63 


0 


5.08 


5.39 


31. 


45 


5.61 


4.99 


3.16 


.24 


5.92 


6.11 



Oi=55°32 . 


0 


1.19 


0 


0 


0 


2.71 


3.64 


33. 


10 


1.65 


0 


0 


0 


3.06 


3.93 


34. 


20 


2.25 


0 


0 


0 


3.56 


4.37 


35. 


30 


3.12 


.79 


0 


0 


4.32 


5.05 


36. 


40 


4.53 


2.60 


0 


0 


5.53 


6.14 


37. 


50 


7.17 


6.15 


3.21 


0 


7.69 


8.00 



39 



1500 



V 2 = 1525 



V L = 



£3 




A/B=l 


A/B=l/2 


A/B=l/5 


A/B=l/10 


A/B=2/l 


A/B=5/l 


38. 


0° 


.69% 


0% 


0% 


. 0% 


2 . 95% 


4 . 34% 


39. ' 


10 


1.19 


0 


0 


0 


3.33 


4.65 


40. 


20 


1.83 


0 


0 


0 


3.88 


5.14 


41. 


30 


2.74 


0 


0 


0 


4.70 


5.90 


42. 


40 


4.22 


.79 


0 


0 


6.02 


7.13 


43. 


50 


6.99 


4.29 


0 


0 


8.38 


9.24 



44. 


0 


0 


0 


0 


0 


3.09 


5.30 


45. 


10 


.11 


0 


0 


0 


3.51 


5.64 


46. 


20 


.79 


0 


0 


0 


4.11 


6.19 


47. 


30 


1.73 


0 


0 


0 


5.01 


7.07 


48. 


40 


3.24 


0 


0 


0 


6.46 


8.49 


49. 


50 


6.07 


.52 


0 


0 


9.07 


10.95 


50. 


55 


12.81 


9.18 


0 


0 


14.71 


15.88 



40 



TABLE 2 

RELATIVE ERROR 
V l = 1525 V 2 - 1500 







A/B=l 


A/B=l/2 


A/B=l/5 


A/B=l/10 


A/B*2/l 


A/B-5/1 


0,=IO O I- 


-40° 


1 . 227c 


07o 


0% 


07. 


1.187. 


2.207. 


1 2. 


-30 


1.27 


.51 


0 


0 


1.64 


1.87 


3. 


-20 


1.31 


.83 


0 


0 


1.55 


1.70 


4. 


-10 


1.40 


1.10 


.20 


0 


1.54 


1.63 


5. 


- 5 


1.46 


1.24 


.58 


0 


1.57 


1.63 


6. 


0 


1.54 


1.39 


.95 


.20 


1.61 


1.65 


7. 


5 


1.64 


1.56 


1.34 


.96 


1.67 


1.70 



0,=2O° 8. 


-10 


1.31 


.83 


0 


0 


1.55 


1.70 


1 9. 


- 5 


1.40 


.99 


0 


0 


1.60 


1.72 


10. 


0 


1.50 


1.17 


.17 


0 


1.66 


1.76 


11. 


5 


1.61 


1.36 


.58 


0 


1.74 


1.82 


12. 


10 


1.76 


1.58 


1.04 


.12 


1.85 


1.90 


13. 


15 


1.93 


1.84 


1.55 


1.06 


1.98 


2.01 



9 =30° 14. 


- 5 


1.37 


.69 


0 


0 


1.71 


1.91 


1 15. 


0 


1.49 


.89 


0 


0 


1.79 


1.97 


16. 


5 


1.63 


1.11 


0 


0 


1.89 


2.05 


17. 


10 


1.79 


1.35 


0 


0 


2.02 


2.15 


18. 


15 


1.99 


1.63 


.54 


0 


2.17 


2.27 


19. 


20 


2.22 


1.96 


1.17 


0 


2.35 


2.43 


20. 


25 


2.50 


2.36 


1.92 


1.19 


2.57 


2.61 



9 =40° 21. 


0 


1.49 


.42 


0 


0 


2.02 


2.33 


1 22. 


10 


1.84 


.93 


0 


0 


2.28 


2.55 


23. 


20 


2.30 


1.60 


0 


0 


2.65 


2.86 


24. 


30 


2.99 


2.55 


1.21 


0 


3.21 


3.34 


25. 


35 


3.46 


3.21 


2.45 


1.16 


3.59 


3.66 



01=50° 26. 


0 


1.37 


0 


0 


0 


2.36 


2.94 


27. 


10 


1.77 


0 


0 


0 


2.66 


3,18 


28. 


20 


2.29 


.68 


0 


0 


3.08 


3.55 


29. 


30 


3.05 


1.72 


0 


0 


3.70 


4.09 


30. 


40 


4.23 


3.35 


.63 


0 


4.66 


4.92 


31. 


45 


5.10 


4.58 


3.00 


.24 


5.35 


5.50 



01=55° 32. 


0 


1.15 


0 


0 


0 


2.55 


3.38 


33. 


10 


1.58 


0 


0 


0 


2.88 


3.64 


34. 


20 


2.14 


0 


0 


0 


3.33 


4.02 


35. 


30 


2.94 


.78 


0 


0 


3.99 


4.61 


36. 


40 


4.18 


2.48 


0 


0 


5.01 


5.51 


37. 


50 


6.34 


5.53 


3.04 


0 


6.73 


6.97 



41 



l x = 1525 



V 2 - 1500 





e 


3 


A/B-l 


A/B=l/2 


A/B=l/5 


A/B=l/10 


A/B=2/l 


A/B-5/1 


0 l =6O° 


38. 


0° 


.66% 


0% 


0% 


07, 


2.73% 


3.95% 




39. 


10 


1.14 


0 


0 


0 


3.08 


4.21 




40. 


20 


1.73 


0 


0 


0 


3.56 


4.63 




41. 


30 


2.58 


0 


0 


0 


4.27 


5.27 




42. 


40 


3.89 


.78 


0 


0 . 


5.38 


6.25 




43. 


50 


6.17 


3.97 


0 


0 


7.23 


7.28 


6 l =65° 


44. 


0 


0 


0 


0 


0 


2.79 


4.67 




45. 


10 


.10 


0 


0 


0 


3.17 


4.95 




46. 


20 


.75 


0 


0 


0 


3.70 


5.40 




47. 


30 


1.63 


0 


0 


0 


4.47 


6.10 




48. 


40 


2.99 


0 


0 


0 


5.66 


7.19 




49. 


50 


5.39 


.51 


0 


0 


7.67 


8.98 




50. 


55 


10.27 


7.81 


0 


0 


11.46 


12.16 



42 



TABLE 3 

RELATIVE ERROR 
* 1525 V 2 - 1475 





03 


A/B=l 


A/B=l/2 


A/B=l/5 


A/B=l/10 


A/B-2/1 


A/B-5/1 


e^io 0 l. 


-40° 


2.43% 


0% 


0% 


0% 


3.65% 


4.36% 


2. 


-30 


2.54 


1.02 


0 


0 


3.28 


3.72 


3. 


-20 


2.63 


1.67 


0 


0 


3.10 


3.39 


4. 


-10 


2.80 


2.21 


.41 


0 


3.09 


3.26 


5. 


- 5 


2.92 


2.49 


1.17 


0 


3.13 


3.26 


6. 


0 


3.08 


2.79 


1.91 


.41 


3.22 


3.31 


7. 


5 


3.27 


3.13 


2.68 


1.93 


3.35 


3.39 



9 1 =20° 8. 


-10 


2.63 


1.67 


0 


0 


3.10 


3.39 


9. 


- 5 


2.80 


2.00 


0 


0 


3.19 


3.43 


10. 


0 


2.99 


2.34 


.34 


0 


3.31 


3.51 


11. 


5 


3.23 


2.72 


1.18 


0 


3.48 


3.63 


12. 


10 


3.51 


3.16 


2.08 


.25 


3.68 


3.79 


13. 


15 


3.85 


3.66 


3.09 


2.14 


3.94 


4.00 



01 = 30 ° 14. 


- 5 


2.75 


1.40 


0 


0 


3.41 


3.80 


15. 


0 


2.98 


1.79 


0 


0 


3.57 


3.92 


16. 


5 


3.26 


2.22 


0 


0 


3.77 


4.07 


17. 


10 


3.58 


2.70 


0 


0 


4.01 


4.26 


18. 


15 


3.95 


3.25 


1.09 


0 


4.30 


4.51 


19. 


20 


4.40 


3.90 


2.35 


0 


4.65 


4.80 


20. 


25 


4.95 


4.67 


3.83 


2.39 


5.09 


5.17 



01=40° 21. 


0 


2.97 


.84 


0 


0 


4.01 


4.62 


22. 


10 


3.65 


1.87 


0 


0 


4.52 


5.03 


23. 


20 


4.56 


3.18 


0 


0 


5.24 


5.64 


24. 


30 


5.88 


5.04 


2.42 


0 


6.30 


6.55 


25. 


35 


6.78 


6.31 


4.85 


2.33 


7.02 


7.16 



0,=5O° 26. 


0 


2.72 


0 


0 


0 


4.64 


5.77 


27. 


10 


3.50 


0 


0 


0 


5.22 


6.23 


28. 


20 


4.53 


1.37 


0 


0 


6.04 


6.92 


29. 


30 


5.98 


3.41 


0 


0 


7.22 


7.94 


30. 


40 


8.21 


6.56 


1.26 


0 


9.02 


9.49 


31. 


45 


9.82 


8.87 


5.89 


.49 


10.29 


10.57 



0j=55° 32. 


0 


2.28 


0 


0 


0 


5.01 


6.58 


33. 


10 


3.13 


0 


0 


0 


5.63 


7.07 


34. 


20 


4.23 


0 


0 


0 


6.49 


7.81 


35. 


30 


5.77 


1.56 


0 


0 


7.75 


8.90 


36. 


40 


8.12 


4.89 


0 


0 


9.66 


10.56 


37. 


50 


12.06 


10.61 


5.96 


0 


12.77 


13.19 



43 



V 1 = 1525 v 2 = 1475 



9 1 =60° 



0 3 A/B=l 



38. 


0°. 


1.307. 


39. 


10 


2.35 


40. 


20 


3.43 


41. 


30 


5.06 


42. 


40 


7.55 


43. 


50 


11.75 



A/B=l/2 A/B=l/5 



0 % 0 % 

0 0 

0 0 

0 0 

1.55 0 

7.72 0 



A/B=l/10 A/B=2/l 

07, 5.32% 

0 5.99 

0 6.92 

0 8.26 

0 10.30 

0 13.63 



0 1 =65° 



44. 


0 


0 


0 


0 


' 0 


5.40 


45. 


10 


.20 


0 


0 


0 


6.13 


46. 


20 


1.48 


0 


0 


0 


7.13 


47. 


30 


3.21 


0 


0 


0 


8.57 


48. 


40 


5.83 


0 


0 


0 


10.77 


49. 


50 


10.31 


1.02 


0 


0 


14.36 


50. 


55 


18.84 


14.64 


0 


0 


20.79 



A/B-5/1 

7.627o 

8.13 

8.91 

10.08 

11.88 

14.73 



8.90 

9.42 

10.26 

11.54 

13.50 

16.62 

21.91 



44 



TABLE 4 

RELATIVE ERROR 
W l = 1500 V 2 * 1550 





£3 


A/B=l 


A/B=l/2 


A/B=l/5 


A/B=l/10 


A/B-2/1 


A/B-5/1 


9]=10° 1. 


-40° 


2.48% 


0% 


0% 


0% 


3.81% 


4.62% 


2. 


-30 


2.57 


1.00 


0 


0 


3.37 


3.85 


3. 


-20 


2.65 


1.65 


0 


0 


3.16 


3.46 


4. 


-10 


2.82 


2.20 


.40 


0 


3.13 


3.32 


5. 


- 5 


2.95 


2.49 


1.14 


0 


3.18 


3.32 


6. 


0 


3.12 


2.81 


1.89 


.39 


3.27 


3.37 


7. 


5 


3.33 


3.17 


2.70 


1.91 


3.41 


3.46 



01=20° 8. 


-10 


2.65 


1.65 


0 


0 


3.16 


3.46 


9. 


- 5 


2.82 


1.98 


0 


0 


3.25 


3.51 


10. 


0 


3.03 


2.34 


.33 


0 


3.38 


3.59 


11. 


5 


3.29 


2.75 


1.15 


0 


3.56 


3.73 


12. 


10 


3.60 


3.21 


2.08 


.24 


3.79 


3.91 


13. 


15 


3.98 


3.77 


3.15 


2.13 


4.08 


4.14 



©,=30° 14. 


- 5 


2.79 


1.38 


0 


0 


3.50 


3.94 


15. 


0 


3.04 


1.79 


0 


0 


3.68 


4.07 


16. 


5 


3.34 


2.23 


0 


0 


3.90 


4.24 


17. 


10 


3.69 


2.74 


0 


0 


4.17 


4.46 


18. 


15 


4.11 


3.33 


1.07 


0 


4.50 


4.74 


19. 


20 


4.62 


4.05 


2.36 


0 


4.91 


5.08 


20. 


25 


5.26 


4.93 


3.97 


2.41 


5.42 


5.52 



CM 

o 

o 

II 

CD 


0 


3.06 


.83 


0 


0 


4.21 


4.91 


22. 


10 


3.80 


1.89 


0 


0 


4.79 


5.38 


23. 


20 


4.83 


3.28 


0 


0 


5.63 


6.11 


24. 


30 


6.40 


5.39 


2.46 


0 


6.91 


7.22 


25. 


35 


7.52 


6.92 


5.17 


2.36 


7.82 


8.00 



01 = 50 ° 26. 


0 


2.85 


0 


0 


0 


5.03 


6.38 


27. 


10 


3.70 


0 


0 


0 


5.71 


6.94 


28. 


20 


4.86 


1.38 


0 


0 


6.68 


7.79 


29. 


30 


6.58 


3.57 


0 


0 


8.15 


9.11 


30. 


40 


9.45 


7.29 


1.26 


0 


10.56 


11.24 


31. 


45 


11.71 


10.34 


6.44 


.48 


12.40 


12.82 



01=55° 32. 


0 


2.41 


0 


0 


0 


5.56 


7.52 


33. 


10 


3.34 


0 


0 


0 


6.29 


8.12 


34. 


20 


4.57 


0 


0 


0 


7.34 


9.06 


35. 


30 


6.38 


1.59 


0 


0 


8.94 


10.53 


36. 


40 


9.39 


5.28 


0 


0 


11.56 


12.91 


37. 


50 


15.22 


12.93 


6.55 


0 


16.40 


17.12 



45 



1500 



V 2 = 1550 



V 



1 “ 





03 


A/B=l 


A/B=l/2 


A/B=l/5 


A/B=l/10 


A/B=2/l 


A/B=5/l 


01=60° 38. 


0° 


1.39% 


0% 


0% 


0% 


6.10% 


9.08% 


39. 


10 


2.42 


0 


0 


0 


6.90 


9.74 


40. 


20 


3.74 


0 


0 


0 


8.06 


10.80 


41. 


30 


5.63 


0 


0 


0 


9.82 


12.47 


42. 


40 


8.76 


1.59 


0 


0 


12.70 


15.20 


43. 


50 


14.84 


8.88 


0 


0 


18.06 


20.08 



01=65° 44. 


0 


0 


0 


0 


0 


10.39 


11.34 


45. 


10 


.22 


0 


0 


0 


11.83 


12.08 


46. 


20 


1.63 


0 


0 


0 


13.90 


13.31 


47. 


30 


3.57 


0 


0 


0 


17.11 


15.29 


48. 


40 


6.72 


0 


0 


0 


22.51 


18.55 


49. 


50 


12.88 


1.04 


0 


0 


33.02 


24.43 


50. 


55 


28.92 


19.98 


0 


0 


60.45 


37.17 



46 



9 1 *10° 



TABLE 5 

RELATIVE ERROR 
V L * 1550 V 2 - 1500 



93 


A/B=l 


A/B*l/2 


A/B»l/5 


A/B=l/10 


A/B*2/l 


A/B-5/1 


1. -40° 


2.39% 


0% 


0% 


0% 


3.59% 


4.29% 


2. -30 


2.50 


1.00 


0 


0 


3.23 


3.66 


3. -20 


2.59 


1.65 


0 


0 


3.05 


3.33 


4. -10 


2.75 


2.18 


.41 


0 


3.04 


3.21 


5. - 5 


2.87 


2.45 


1.15 


0 


3.08 


3.21 


6. 0 


3.03 


2.74 


1.88 


.40 


3.17 


3.25 


7. 5 


3.22 


3.07 


2.64 


1.90 


3.29 


3.34 



O^O 0 8. 


-10 


2.59 


1.65 


0 


0 


3.05 


3.33 


9. 


- 5 


2.75 


1.96 


0 


0 


3.14 


3.37 


10. 


0 


2.94 


2.30 


.33 


0 


3.26 


3.45 


11. 


5 


3.17 


2.68 


1.16 


0 


3.42 


3.57 


12. 


10 


3.45 


3.10 


2.05 


.24 


3.62 


3.73 


13. 


15 


3.79 


3.60 


3.04 


2.10 


3.88 


3.93 



9 i =30° 14. 


- 5 


2.70 


1.37 


0 


0 


3.35 


3.74 


15. 


0 


2.94 


1.76 


0 


0 


3.51 


3.85 


16. 


5 


3.20 


2.19 


0 


0 


3.71 


4.00 


17. 


10 


3.52 


2.66 


0 


0 


3.94 


4.20 


18. 


15 


3.89 


3.20 


1.07 


0 


4.23 


4.43 


19. 


20 


4.33 


3.83 


2.31 


0 


4.58 


4.73 


20. 


25 


4.87 


4.60 


3.77 


2.35 


5.01 


5.09 



9, =40° 21. 


0 


2.93 


.83 


0 


0 


3.94 


4.55 


22. 


10 


3.59 


1.84 


0 


0 


4.44 


4.95 


23. 


20 


4.49 


3.13 


0 


0 


5.16 


5.55 


24. 


30 


5.79 


4.96 


2.38 


0 


6.20 


6.45 


25. 


35 


6.68 


6.21 


4.78 


2.29 


6.91 


7.05 



9!=50° 26. 


0 


2.68 


0 


0 


0 


4.57 


5.68 


27. 


10 


3.45 


0 


0 


0 


5.14 


6. 14 


28. 


20 


4.46 


1.35 


0 


0 


5.94 


6.82 


29. 


30 


5.89 


3.36 


0 


0 


7.11 


7.82 


30. 


40 


8.09 


6.46 


1.24 


0 


8.88 


9.35 


31. 


45 


9.67 


8.73 


5.80 


.48 


10.14 


10.41 



9 t =55° 32. 


0 


2.24 


0 


0 


0 


4.93 


6.48 


33. 


10 


3.08 


0 


0 


0 


5.54 


6.96 


34. 


20 


4.16 


0 


0 


0 


6.39 


7.69 


35. 


30 


5.68 


1.53 


0 


0 


7.63 


8.77 


36. 


40 


8.00 


4.81 


0 


0 


9.51 


10.40 


37. 


50 


11.89 


10.45 


5.87 


0 


12.59 


13.00 
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V]_ = 1550 



1500 



V 2 = 





03 


A/B=l 


A/B=l/2 


01=60° 38. 


0° 


1.28% 


0% 


39. 


10 


2.21 


0 


40. 


20 


3.37 


0 


41. 


30 


4.98 


0 


42. 


40 


7.44 


1.52 


43. 


50 


11.57 


7.60 



A/B=l/5 


A/B=l/10 


A/B=2/l 


A/B=5/l 


0% 


0% 


5.24£ 


7.50% 


0 


0 


5.90 


8.00 


0 


0 


6.81 


8.77 


0 


0 


8.13 


9.33 


0 


0 


10.15 


11.71 


0 


0 


13.44 


14.52 



0 1 =65° 44. 


0 


0 


0 


0 


0 


5.32 


8.77 


45. 


10 


.20 


0 


0 


0 


6.03 


9.29 


46. 


20 


1.46 


0 


0 


0 


7.02 


10.11 


47. 


30 


3.16 


0 


0 


0 


8.44 


11.37 


48. 


40 


5.74 


0 


0 


0 


10.61 


13.31 


49. 


50 


10.16 


1.01 


0 


0 


14.16 


16.39 


50. 


55 


18.58 


14.53 


0 


0 


20.51 


21.63 
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TABLE 6 

ABSOLUTE ERROR 



V 1 = 1500 V 


2 * 1525 






Vi » 


1525 V 2 


- 1500 




®3 


A/B=l 


A/B=5/l 


A/B=l/5 


A/B=l 


A/B=5/l 


A/B-l/5 


0!=1O° 1. 


-40° 


.03% 


.23% 


0% 


.03% 


.23% 


0% 


2. 


-30 


.03 


.26 


0 


.04 


.26 


0 


3. 


-20 


.05 


.32 


0 


.05 


.32 


0 


4. 


-10 


.08 


.53 


0 


.08 


.47 


.03 


5. 


- 5 


.11 


.62 


.02 


.11 


.62 


.09 


6. 


0 


.17 


.94 


.08 


.18 


.96 


.21 


7. 


5 


.37 


“ 


.22 


.32 


— 


.61 


01=30° 8. 


- 5 


.04 


.30 


0 


.04 


.30 


0 


9. 


0 


.05 


.36 


0 


.06 


.36 


0 


10. 


5 


.07 


.44 


0 


.07 


.45 


0 


11. 


10 


.10 


.58 


0 


.09 


.58 


0 


12. 


15 


.14 


.80 


.08 


.15 


.81 


.07 


13. 


20 


.23 


1.28 


.24 


.24 


1.28 


.24 


14. 


25 


.52 


2.29 


.79 


.51 


2.73 


.80 


01=50° 15. 


0 


.03 


.30 


0 


.03 


.30 


0 


16. 


10 


.04 


.39 


0 


.04 


.39 


0 


17. 


20 


.07 


.56 


0 


.07 


.55 


0 


18. 


30 


.13 


.93 


0 


.12 


.92 


0 


19. 


40 


.35 


2.07 


.10 


.32 


2.09 


.10 


20. 


45 


.82 


4.46 


.95 


.83 


4.50 


.96 



49 



BIBLIOGRAPHY 



i # Bush, V. and Conant, J. B. The Application of Oceanography 
to Subsurface Warfare (Military Oceanography) . National 
Defense Research Committee, 1946. 



2. 


Encyclopedia Americana, v. 13, 


20, 1964. 




3. 


Kinsler, L. E. and Frey, A. R. 
2nd ed. John Wiley and Sons, 


Fundamentals of 
1962. 


Acoustics , 


4. 


Thompson, B. and Anderson, R. 


Subjective Layer 


Depth 



Analysis Model . U. S. Naval Oceanographic Office, 1964. 




50 



n n n n n n 



Appendix I 



PROGRAM THESIS3 

ASSUME TWO BUOYS A DISTANCE OF 2A APART. ORIGIN (0*0) IS INBETWEEN 
ACTUAL TARGET IS AT(X3»Y3)» APPARENT TARGET ATIX4.Y4). ASSUME 2 
DIFFERENT VELOCITY GRADIENTS VI AND V2 WHERE V2 IS A DISTANCE B FROM 
THE BUOY AXIS. 

READ 1000 * K » VI » V2 
1000 FORMAT ( I5*F10.3*F10.3) 

HOLD EITHER THETA 1 OR THETA 3 CONSTANT. K IS NUMBER OF VARYING 
THETA OR A OR B TO BE READ IN. 

J = 0 
L=K+1 

7 J=J+1 

I F ( J.EO.L) 1*2 

2 READ 2000»A,3*THETA1»THETA3 
2000 FORMAT (F6.3»F6.3*F10.6»F10.6 ) 

PRINT 5000 * VI » V2 
5000 FORMAT ( F10.3»F10.3) 

PRINT 200*A.B*THETA1»THETA3 

200 FORMAT ( 1H1 19X > 5H A = F6.3*10X,5H B = F6.3*I0X*10H THETA1 « F10.6* 
*10X * 10H THETA3 = F10.6//) 

I F ( VI .LT. V2 ) 3 *4 

3 THETAC=ASINF(V1/V2 ) 

PRINT 300 »THET AC 

300 FORMAT ( 20X ♦ 10H THETAC = F10.6//) 

I F ( THETA1 .LT. THETAC ) 5 * 6 

5 CONTINUE ! 

4 TT1=TANF{ THETA1 ) 

TT3=TANF(THETA3) 

X4=2*A/(TT1-TT3 ) 

Y4=A*(TT1+TT3)/(TT1-TT3) 

PRINT 400*TT1»TT3 

400 FORMAT ( 20X ♦ 7H TT1 = F9.6//20X*7H TT3 = F9.6//) 

PRINT 500 »X4»Y4 

500 FORMAT ( 20X ♦ 22H APPARENT POSITION = <F5.2*2H *F5.2»2H )//) 
Z1=ASINF(V2/V1*SINF( THETA1 ) ) 

Z3=ASINF(V2/V1*SINF( THETA3 ) ) 

T 1 = TANF { Z 1 ) 

T3=TANF(Z3) 

PRINT 600, Z1 ,Z3»T1 »T3 

600 FORMAT ( 20X , 6H Z1 = F9.6,10X,6H Z3 = F9.6*10X*6H T1 ■ F9.6*10X* 

*6H T3 = F9.6//) 

X3= <8*(TT3-TT1+T1-T3)+2*A )/ (T1-T3) 

Y3=(B*( T3*TT1-T1*TT3 )-A*( T1+T3) )/( T3-T1 ) 

IF ( X3 .LE. B ) 8 * 9 

8 X3=X4 
Y3=Y4 

9 PRINT 700 »X3 » Y3 

700 FORMAT ( 20X » 18H REAL POSITION .= (F5.2,2H ,F5.2»2H )//) 

E= ( SORTF ( (X3-X4)**2+(Y3-Y4)**2) ) * 1 00 .0/ SORTF ( X3**2+Y3**2 ) 

PRINT 800, E 

800 FORMAT ( 20X, 16H ERROR VECTOR <= F5.2, 8H PERCENT) 

GO TO 7 

6 PRINT 900 

900 FORMAT ( 20X > 53H ANGLE GREATER THAN CRITICAL ANGLE. TOTAL REFLECTION 
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*■ • ) 

GO TO 7 
1 CONTINUE 
END 

END 

FINIS 

XECUTE » 
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